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Abstract 

Exact exchange (EXX) Kohn-Sham calculations within an all-electron full-potential method are 
performed on a range of semiconductors and insulators (Ge, GaAs, CdS, Si, ZnS, C, BN, Ne, Ar, Kr 
and Xe). We find that the band-gaps are not as close to experiment as those obtained from previous 
pseudopotential EXX calculations. Full-potential band-gaps are also not significantly better for 
sp semiconductors than for insulators, as had been found for pseudopotentials. The locations of 
d-band states, determined using the full-potential EXX method, are in excellent agreement with 
experiment, irrespective of whether these states are core, semi-core or valence. We conclude that 
the inclusion of the core- valence interaction is necessary for accurate determination of EXX Kohn- 
Sham band structures, indicating a possible deficiency in pseudopotential calculations. 

PACS numbers: 71.15.Mb,71.10.-w,71.22.+i 



Exact treatment of exchange within the Kohn-Sham (KS) formulation of density func- 
tional theory (DFT) has been one of the most interesting developments in the last few years 
The KS band-gaps (E^ s ) for sp semiconductors, determined using exact 
exchange (EXX) and the pseudopotential (PP) method, are found to be very close to ex- 
perimental gaps 0,0], even in the absence of any correlation. This is surprising because the 
relation between E^ s and the fundamental band-gap E g of solids is given by E g = E^ s + A xc 
where, A xc is the discontinuity in the exchange-correlation potential. This discontinuity is 
lacking in the local density approximation (LDA), resulting in the equality of the KS and 
fundamental band-gap 5} which leads to the famous band-gap problem. But for the EXX 
potential even in the absence of exact correlation, the discontinuity does not vanish and 
hence one cannot expect that the EXX KS gap alone can be identified directly with ex- 
perimental values. In fact, Stadele et al. [3J determined the discontinuities in the exchange 
potential explicitly for Si, Ge and GaAs, and found them to be about 3-5 times the band-gap. 
Furthermore, Magyar et al. j^] noted that the EXX KS gap for wide band-gap noble gas 
solids is not as close to the experimental gap (either optical or fundamental) as it is for sp 
semiconductors. One of the major conclusions of this work was that the EXX method does 
not provide a KS band structure that agrees equally well with experiment for semiconductors 
and insulators. All the above work has been done using the PP approach. This approach 
neglects both the exact core-valence exchange interaction as well as the relaxation of the 
core states in the solid-state environment. There do exist all-electron EXX calculations of 
the band-gaps in C, Si and Ge 4]. In this work, the linear combination of muffin-tin orbitals 
(LMTO) within the atomic sphere approximation (ASA) used indicates that the band-gaps 
are an improvement over LDA but do not agree with experiment as well as PP-EXX. These 
calculations however suffer from the shape approximation of the potential which ignores 
non-spherical contributions. In our recent article, the importance of this non-spherical part 
of the self-consistent EXX potential was demostratedj^]. 

At this point it is relevant to note that another method used to improve upon the single 
particle KS band-gaps is the GW approximation. Here a partially or fully self-consistent so- 
lution of the Dyson equation is used to determine the self-energy which provides a correction 
to the KS-LDA band-gap. It has been inferred from GW calculations performed using the 
PP approach that full self-consistency, even though it is essential for charge conservation, 
worsens the agreement of quasi-particle band-gaps 
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recently it was noted by Ku and Eguiluz Q| that the success of the partially self-consistent 
GW approach is due to some kind of error cancellation which is an artifact of the PP and 
also that the deep-core states have an important effect on the band-gaps. Such an effect 
may be of similar importance for EXX calculations and requires thorough investigation. 

Another important contributory factor to the gap error is the cation <i-band position in 
sy semiconductors, especially when these levels are relatively shallow (e.g. GaAs, ZnS, CdS) 



12j. In these cases, use of the LDA leads to a misplacement of the <i-bands, although this 
can be improved upon by the GW approximation or the self-interaction correction (SIC) 
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Previous PP-EXX results indicate that even though the EXX 



potential is self-interaction free and the band-gaps are an improvement over those of LDA, 
the (i-band positions are still not correct [l5l |. 

In this Letter, we present results of all-electron full-potential EXX calculations applied to 
representative semiconductors and insulators. This is to ascertain the effect of core-valence 
interaction and asymmetry in the potential on band-gaps and <i-band positions, and also to 
find the reason for the inconsistency in EXX band-gaps of semiconductors and insulators. 

We have implemented EXX potential using the all-electron full-potential linearized 
augmented-plane wave (FP-LAPW) method Q within the EXCITING code Q- The start- 
ing point for these calculations is the exchange energy 
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where tp^ is the zth Kohn-Sham orbital of fc-point k, is the /c-point weight, n(r) = 
Sik C u, k0*k( r )0«k( r ) is the density, and i and j run over both core and valence states, in 
contrast to PPs where the indices run only over the valence states. Proceeding in the same 
manner as Gorling et al. , the functional derivative chain rule is employed to obtain 

the exchange potential 
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where are the Kohn-Sham eigenvalues and 
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In order to obtain the functional derivative Sv s (r') / 5n(r) , the linear- response operator Xi 
defined by the relation 



must be inverted. This is not possible directly as x has a zero eigenvalue corresponding to a 
constant eigenfunction. However, in a basis which excludes such functions the inversion of 
X can be performed. Application of elementary perturbation theory to n leads to an explicit 
expression for the response 



Implementations of EXX for pseudopotentials have used the plane wave basis for expanding 
X- This is not feasible in a FP-LAPW approach because of the strongly varying wave 
functions inside the muffin-tin. Instead, we propose a new basis as follows. The overlap 
densities p a (r) = 0* k (r)0, k (r) and their complex conjugates are used as a spanning set, where 
a = (ik, jk) with i and j labelling occupied and unoccupied states, respectively. These are 
the states of some representative /c-point which is chosen in advance. This spanning set is 
then reduced to a basis set as follows. First the overlap matrix of the spanning set elements 



is obtained and diagonalized. Note that O, which is positive semidefinite, has nonnegative 
eigenvalues. Next, all the eigenvectors of O corresponding to eigenvalues smaller than a 
certain tolerance, e, are discarded. Finally, denoting the remaining eigenvectors as v@, 
where (3 labels the vector and 7 its coordinate, we seek a transformation matrix C such that 
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then 




Straight-forward algebra shows that C should satisfy 



CC ] = (v^Ov) 
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where v is the matrix of eigenvectors in column-wise form. The matrix C is obtained by 
performing a Cholesky decomposition on the right hand side of Eq. |3] By virtue of their 



basis for the expansion and inversion of Xi as we h as for the term in square brackets in Eq. 
[0 We should also point out that this basis may be useful for FP-LAPW time-dependent 
DFT response and GW methods. 

Special attention should now be drawn to the calculation of the non-local matrix elements 
(Eq. |2J) which may be determined by a well-established method for solving Poisson's equation 
in a FP-LAPW environment Q]. The differences here are that the densities are now complex 
and there is a long-range term arising when q = k — k' is close to zero. This is treated by 
considering the so-called pseudo-charge density Q], which is chosen to be sufficiently smooth 
within the muffin-tins so that it may be expanded in terms of plane waves, and yet has the 
same multipole expansion as the real density. If we restrict k and k' such that q is in the 
first Brillouin zone (BZ), then the long-range (LR) contribution to the matrix elements from 
the pseudo-charge is 



where pu(q) and pzj(q) are the pseudo-charge densities in reciprocal space. This sum suffers 
from poor convergence with respect to the number of g-points. We therefore approximate it 
by an integral over a sphere of volume equivalent to that of the BZ. The final expression is 



All the calculations in the present work are performed on a mesh of 27 special Appoints 
in the irreducible BZ. As the response and its inverse are very sensitive to the number of 
empty states, we have taken special care to use a sufficient number to converge the band- 
gaps to within 0.01 eV. We find that almost all the materials studied here have converged 
d-band positions and band-gaps with around 25 empty states. All the results presented are 
generated using 30 empty states. One final practical note is that the eigenvalue cut-off, e, 
discussed above, is taken to be 10 -5 times the largest eigenvalue of O. For example for Si 
this results in 73 elements in the full basis set. 

In Fig. 1 we present the difference between the experimental fundamental band-gap and 
the KS band-gap for various semiconductors and insulators. The EXX potential is self- 



construction and the fact that J dr p a (r) = 0, the set of functions {p a } form an optimal 
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FIG. 1: Difference between theoretical KS gaps and experimental fundamental gaps in eV. Full- 
potential exact exchange (FP-EXX), FP-EXX without core- valence interaction (FP-EXX-ncv), 
FP-EXX with only spherical exchange (FP-EXX-sph) and LDA (FP-LDA) results are from this 
work. The PP-EXX results are from Refs. 5J, d, 15. The LMTO-ASA-EXX results are from Ref. 



and the experimental data is taken from Refs. [IS, 
is -10.25 eV, and is off the scale of the plot. 
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2Cl Note that LDA gap difference for Ne 



interaction free and hence an enhancement of the KS band-gap over that of the LDA is 
expected. This is indeed the case here. We find that for semiconductors the FP-EXX band- 
gaps are overestimated with respect to experiments by up to 60%. The band-gaps in all 
the insulators except Xe are underestimated with respect to the experimental fundamental 
gap up to about 24%. For Kr and Xe the KS band-gap is overestimated with respect to 
the experimental optical gap by 3.8% and 24.6% respectively, while for other insulators it is 
underestimated up to 6% . 
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Unlike for PP-EXX, there is not the excellent agreement of the KS gaps of sp semicon- 
ductors with experiment. In contradistinction to the PP-EXX results we find that the EXX 
is not inconsistent in its treatment of semiconductors and insulators, in fact agreement with 
experiment of the KS gaps for insulators is not any worse than for the semiconductors. In 
order to determine the reason for this behaviour we performed ground state calculations 
for these materials by ignoring the core-valence interaction term in the EXX potential. We 
emphasise here that this term was removed only from the exchange potential and not from 
the rest of the potential. Even this has substantial effect on the KS band-gaps. The results 
are marked as FP-EXX-ncv in Fig. 1. Compared to FP-EXX, the FP-EXX-ncv band-gaps 
for all the semiconductors are closer to the PP-EXX and the experimental values. On the 
other hand, the agreement of the KS band-gaps with experiments for insulators now wors- 
ens. This trend is concomitant with the findings of the PP-EXX calculations 0, 0] where 
the performance of the EXX potential in determining the KS gaps is excellent for the semi- 
conductors and lacking for the insulators. Previous all electron EXX calculations performed 



within the ASA show a similar oyorestnuatiou of the gaps for semiconductors Q. Ail this 
suggests that the lack of core-valence interaction is responsible for the spurious agreement 
of the PP-EXX KS band-gaps with the experiments for sp semiconductors. 

In our work on magnetic metals [7| we found that inclusion of the non-spherical contribu- 
tions to the exchange potential is crucial to obtain the correct ground-state. The results for 
the KS band-gaps obtained using only the spherical contribution to the EXX potential are 
marked as FP-EXX-sph in Fig. 1. The effect of the shape approximation to the potential 
in all the materials other than noble gas solids is substantial. Which is expected since the 
noble gas solids are composed of almost independent atoms. 
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Another property which is known to lead to errors in the band-gaps 
LDA performs badly is the position of the semi-core (i-states. The LDA under-binds leading 
to displacement of <i-bands much above the experimental value. In the past, a quasi-particle 
bandstructure was determined u sing the self energy calculated with the GW approximation 
to correct the <i-band position |l3l| . Also the SIC to the potential was found to lead to 
good agreement between the theoretical and experimental position of these (i-bands. Since 
the EXX is self-interaction free in nature, it should improve the <i-band positions. Counter 
intutively Rinke et al. found that even though the PP-EXX calculations result in a 
downward motion of the <i-bands with respect to LDA, still GW calculations are needed to 
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FIG. 2: Difference between theoretical (i-band eigenvalues and experimental data in eV. Full- 
potential exact exchange (FP-EXX), FP-EXX without core- valence interaction (FP-EXX-ncv) and 
LDA (FP-LDA) results are from this work. The results obtained using the GW calculation on 
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the self-interaction corrected 



top of the PP-EXX ground-state (GW-PP-EXX) are from Ref. 
potential within PP approach (PP-EXX) are from Ref. 

LMTO method (GW-FPLMTO) and Hartree-Fock FP-LMTO (HF-FPLMTO) are from Ref. 
The GW calculation on top of the LMTO-ASA-LDA ground-state (GW-LMTO-ASA) are from 
Ref. 



14L GW results calculated using a FP- 



fully correct the positions. The FP-EXX potential, in which the core and valence states are 
treated on the same footing, is truly self-interaction free in nature and it would be interesting 
to see how this effects these semi-core/core-like rf-bands. 

In the present work, we compare FP-EXX <i-band eigenvalues with those of previous 
PP-EXX, experiment and various GW and SIC calculations in Fig. 2. The compounds are 
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chosen so that <i-bands are in various energy regimes to get an overall picture: core-like 
(Ge: experimentally -29.6 eV from the top of the valence band), semi-core (GaAs: -18.8 eV, 
InP: -16.8 eV) and valence (ZnS: -8.7 eV, CdS: -9.2 eV). As can be seen, the eigenvalues of 
d-states obtained using the FP-EXX method are in excellent agreement with experiments. 
In fact, in several cases it is better than that of the FP-GW results. 

Here again the agreement between the FP-EXX-ncv and PP-EXX values of Rinke et al. 
is quite good. This indicates that the inclusion of core-valence interaction is crucial for 
determination of the correct <i-band eigenvalues. 

The main aim of our paper is to establish the effects of the inclusion of core-valence 
interaction in full-potential EXX calculations. Prior to this investigation, it appeared that 
EXX might be the answer to the band-gap problem at least for sp semiconductors. We 
find that this is not the case if core-valence interaction is included in the calculation of the 
EXX potential. We also find that EXX is consistent in its treatment of semiconductors and 
insulators. Our findings indicate a possible shortcoming in the PP-EXX method arising from 
the lack of explicit core- valence interactions in the PP-EXX potential, as mentioned earlier by 
Ku and Eguiluz d|. Furthermore, we show that the <i-band eigenvalues determined using 
the FP-EXX method are in excellent agreement with experiment irrespective of whether 
these states are in the core, semi-core or valence. In most cases the agreement is better than 
that achieved by the full-potential GW calculations. 
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